Entanglement of mechanical oscillators coupled to a non-equilibrium environment 
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Recent experiments aim at cooling nanomechanical resonators to the ground state by coupling 
them to non-equilibrium environments in order to observe quantum effects such as entanglement. 
This raises the general question of how such environments affect entanglement. Here we show that 
there is an optimal dissipation strength for which the entanglement between two coupled oscillators 
is maximized. Our results are established with the help of a general framework of exact quantum 
Langevin equations valid for arbitrary bath spectra, in and out of equilibrium. We point out why 
the commonly employed Lindblad approach fails to give even a qualitatively correct picture. 
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1. INTRODUCTION 

Entanglement [1 constitutes a cornerstone of quantum 
mechanics and is a major subject of present-day research 
Whether it persists and can be observed in systems 
comprising macroscopic bodies has been a hotly debated 
topic since the early days of quantum mechanics. The 
ground state of two interacting quantum systems will 
generically be entangled. Thus, one could naively ex- 
pect that it is sufficient to simply cool two interacting, 
macroscopic bodies to their ground states and thereby 
prepare an entangled state. However, when coupling to 
a dissipative bath - as is of course necessary for cooling - 
entanglement may be destroyed, as explored in a number 
of works, e.g., [3 . A slate of recent experiments has now 
brought a new aspect into focus: A non- equilibrium envi- 
ronment, consisting of either a driven optical cavity [4 , a 
superconducting microwave resonator [5] or a supercon- 
ducting single electron transistor [6^, can be employed 
to cool the motion of mechanical resonators down to the 
ground state. The advances in this field may ultimately 
enable tests of quantum mechanics in an entirely new 
regime \T and to observe entanglement of massive ob- 
jects [HI [9]. Still it remains to resolve the issue of how 
the dissipative coupling to the non-equilibrium bath af- 
fects entanglement. 

In the present work, we demonstrate a non-monotonic 
dependence of entanglement between two oscillators on 
the coupling strength to the non-equilibrium environ- 
ment and show that there is an optimal value for the 
coupling to the bath. Below this value, entanglement is 
diminished by thermal fluctuations, and above this value, 
it is lost through dissipation. The striking behavior found 
here is missed entirely by the commonly employed Lind- 
blad approach to dissipative dynamics. 



In order to obtain an exact description, we develop 
a general framework based on quantum Langevin equa- 
tions, which allows us to analyze the entanglement be- 
tween harmonic oscillators in the presence of coupling to 
a linear bath of arbitrary spectral density. First, we ex- 
ploit this scheme to show that even in equilibrium there 
are effects likely to be missed by simpler approaches. 
For example, the minimum coupling strength needed for 
entanglement depends logarithmically on the cutoff fre- 
quency for the most important case of an Ohmic bath 
spectrum. For the case of a non-equilibrium bath, we 
illustrate the generic behavior in a concrete example of 
two mechanical resonators inside an optical cavity, being 
cooled by the optomechanical interaction with the light 
field circulating in the cavity. 




Figure 1. The system consists of two harmonic oscillators {A 
and B) that are coupled via a harmonic force with spring con- 
stant k and are subject to fluctuating quantum forces (Fa/b) 
due to their coupling to the environment. 



2. MODEL 

We consider two coupled oscillators with masses ttia^b 
and frequencies ^a,b^ cf. Fig. [T] In terms of their po- 
sitions and momenta, Qa/b and Pa/b-, the Hamiltonian 
reads ^sys = Ea=A,B ^«^a<7a/2 ^ pl/2ma + k{qA - 
to)^/2, with a coupling spring constant k. Moreover, 
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4. EQUILIBRIUM BATH 



we assume the oscillators to be subject to fluctuating 
quantum forces Fq,, which are possibly correlated, and 
which derive from a bath of harmonic oscillators, with 
^sys-bath = Y^a^ocFa- They wih be characterized by 
their spectra as specifled below. 

If the state of the environment is Gaussian, the oscil- 
lators also end up in a Gaussian state, which is fully de- 
scribed by the covariance matrix = tr(p{i?i, i?j}/2). 

Here R = {paAa-,PbAb)^ = in steady state, 
and p is the system's density matrix. As a measure of 
the entanglement between the oscillators, the logarithmic 
negativity p!QHT2] is calculated as En{p) = Xli=i 2 fi^'i)^ 
where /(c) = — log2(2c) for c < 0.5 and /(c) = oth- 
erwise, and where Ci^2 are the symplectic eigenvalues of 
the partially transposed covariance matrix 7-^^ [12 . 

For later use, and in order to flx the notation it will 
be convenient to consider first the simple example of two 
identical oscillators {rriA/B = ^7 ^a/b = ^) at thermal 
equilibrium and assume the coupling to the environment 
to be negligible. The system can be decoupled by intro- 
ducing the normal- mode coordinates fi± = {qa :^qB)/V^ 
and momenta 7r± = {pA :^Pb)/V^ corresponding to the 
center-of-mass motion (1]+ = Q) and the relative motion 
at frequency = ^yQ'^ -\- 2k/m ^ Vt^ 2G. Here we de- 
fined the coupling rate G = k/2mQ^ to be used in place 
of k. In the following, for simplicity, we assume attrac- 
tive interaction, G oc k > 0. The symplectic eigenvalues 
have the simple form 



Cl,2 



(1) 



and the variances are given by (7)^) = {2n± + l)/2mQ± 

and (tt^) = m(^±(2n± + 1)/2, where n± = (e^±/^-l)"^ 
is the thermal occupation number (we set ks = 1 and 
h = 1). Entanglement is obtained when the product 
(tt^) in Eq. ([l]) becomes smaller than 1/4, which 
requires a coupling rate G/Cl ^ 2nth. In this expression 
only terms up to first order in G/Q have been consid- 
ered and we have set nth = ^+ ^ For a given 
coupling rate G, the logarithmic negativity decreases 
linearly as a function of the thermal occupation nth- 
En{p) - {2G/n - 4nth)/ln2 (black curve in Fig. [§). 
Thus, as is well known, thermal fiuctuations will reduce 
and eventually destroy entanglement. 



eliminating the bath degrees of freedom: 

rriaQait) + rria^lqait) + k{qa{t) - qa{t)) = 
t 

Fa{t)+ I X^p{t-t')qp(t')dt', 



(2) 



(5=A,B 



where a = A/B and a = B/A. Fj^/b denotes 
stationary quantum noise forces acting on the oscilla- 
tors (with (^fJ^ = 0). The response functions that 
take into account the memory effect of the baths 
are given by x^/sW = "^^W ( [^a W, ^/sW] )• 
Solving Eq. Q in Fourier space yields posi- 
tion correlators {qaQ(3)uj = / dte^^'^ {qa{t)qf3{{))) = 
Ec.,^G{A,B} XaaHXf3^{-^){hF^)uj- Here 
{FaF^)o. = Jdte'^'{F^{t)F^{0)) and Xa^{oj) are 
elements of a matrix whose inverse is given by 
U"'H)c.c. = ^«(^' - UJ^) ^ k - xLM and 
(x"^(cj))^^ = -k - Xafdi^) a ^ p. Momentum 
correlators follow from {paPi3)^ — vfio.m^ijp' (^a^/?)^, and 
{Paqis)^ = —irriaUJa {qaq^)^' Finally, equal-time correla- 
tors are obtained by integration, {qaq(3) = / {Qaq^d)^- 
The solution of Eq. Q thus provides the full covariance 
matrix 7 in terms of frequency integrals over arbitrary 
bath spectra, and from it the logarithmic negativity En 
for two coupled dissipative oscillators. Note that we did 
not assume equilibrium, i.e., the fiuctuation-dissipation 
relation between Xa/s does not necessarily 

hold. 

For simplicity, we will from now on restrict our explicit 
calculations to the symmetric case of two identical oscil- 
lators that couple equally strongly to independent baths, 
such that (FaFf^)^ = 6^(3 (FF)^ and Xa(3 = ^a(3X^ - 
The system can then, as before, be decomposed into the 
center of mass mode (7)+, 7r+) and the relative mode 
(/}_, 7r_), which become independent dissipative oscil- 
lators. We find {fj±fj±)^ = (^^^^ |x±('^)l^7 where 

X±{^) — {vfi{S^± — bp') — x^ (bp)) ^ . After frequency 
integration, Eq. ([T]) thus directly yields the logarithmic 
negativity. 



4. EQUILIBRIUM BATH 



3. EXACT SOLUTION 

Returning to the full model, an exact description of the 
dissipative dynamics is provided by quantum Langevin 
equations [13j for the Heisenberg operators, obtained by 



First, we illustrate the general scheme for the case 
of equilibrium baths, picking the important example 
of an Ohmic bath spectrum: {FaFP)^ = (FF)^ = 
mrmOo{coth{bj/2T) + 1)/(1 + a;V^?)- Here T denotes 
the temperature, F^ the damping rate, and Uc the cut- 
off frequency. For F^ < Q and bJc ^ the position 
and momentum variances of an oscillator coupled to this 
bath are given analytically in [14^. Here we only display 
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thermal occupation nth 



Figure 2. Entanglement between two identical harmonic os- 
cillators as a function of the thermal occupation number nth 
in the absence of dissipation (black curve) and under the in- 
fluence of an equilibrium Ohmic bath (red curve), which leads 
to a temperature-independent reduction of entanglement by 
Fm ( In ^ — l)/ (ttQ In 2). The additional reduction with tem- 
perature scales as at low T, as shown in the inset (dashed 
curve). G 0.2 Q, Fm = O.IQ, and ujc 10 

the expansion to first order in F^/Q at T = 0: 

2<*|>/mO,.l + -r^(21n^-l). (3) 

As illustrated in Fig. |2j entanglement between the os- 
cillators is suppressed due to their coupling to the bath. 
The high-frequency bath modes cause momentum fluc- 
tuations that depend logarithmically on the cutoff fre- 
quency, cf. Eq. ([3|. Thus, even at zero temperature, 
the coupling to the environment reduces the logarithmic 
negativity by Fm{^^^ — l)/('7rr21n2) [as follows from 
Eqs. ([T]) and ([3|], and eventually destroys the entangle- 
ment completely. Entanglement persists {Ejy > 0) only 
if the coupling rate exceeds a threshold value of 

Gn.i„ =— (ln--l). (4) 

As a distinctive feature, the minimal coupling rate de- 
pends logarithmically on the cutoff frequency. It indi- 
cates that any approach that disregards the influence of 
high-frequency fluctuations has to fail, as discussed for 
the example of the Lindblad approach below. Our gen- 
eral formula also allows to obtain the full temperature- 
dependence, cf. Fig. [2] 

5. NON-EQUILIBRIUM BATH 

Tunable non-equilibrium quantum fluctuations are 
now relevant in many contexts, and may be used, e.g., 
to cool systems below the bulk temperature. A paradig- 
matic example is the photon shot noise coupled to me- 
chanical resonators in optomechanical setups ^ ^15 (the 



following results also apply to analogous electromechan- 
ical systems [SI [6]). We treat the conceptually clearest 
case where two nanomechanical membranes are placed 
inside a laser-driven cavity and two independent light 
forces F^"^ act on the mechanical normal modes fj± 
leading to optomechanical cooling [161 Hi]- This may 
be realized in a setup with two cavity modes, where 

^sys-cav = (^/-^m) ((«+ + «+)^+ + + ^- ) , cf. Fig. 

|3] Here a± are the annihilation operators of the cavity 
modes, £rn = l/V'^rafl is the mechanical ground-state 
width, and g is the oscillator-cavity coupling rate that 
scales linearly with the laser amplitude (see [18l [19] for a 
derivation of this type of coupling) . The mechanical cou- 
pling k between the oscillators (here assumed as given) 
can itself be implemented via other, strongly driven far- 
detuned cavity modes [H |18] . Other possible setups in- 
clude cold-atom or hybrid atom-membrane systems [18]. 




Figure 3. Two coupled mechanical oscillators, represented 
by membranes, under the influence of non-equilibrium pho- 
ton shot noise baths because of their coupling to two modes 
(a±) of an optical cavity. The placement of the membranes 
allows the normal modes of the coupled system to be cooled 
by independent noise forces F±^^ . 

Elimination of the cavity degrees of freedom generates 

/ ^ ^ \ cav 

cavity noise spectra [T6l (F±F±) = (^/^^)^/^[(co' + 

A±)^ + where hi is the decay rate of the cav- 

ity photons and A± the detuning of the correspond- 
ing input lasers with respect to the flrst (second) cav- 
ity mode. A spectrum of this kind induces an op- 

tomechanical cooling rate of ropt,± = ^myY±^±/ ^ ~ 
(^F±F±^ y In the optimal cooling regime, for A± = 

—Q± and <C 1^, we have Fopti ^ Fopt = 4^^/^- 
In this regime, the minimum possible phonon number 
due to optical cooling, deflned by (nopt± + l)/^opt± = 

/ ^ ^ \ cav / ^ ^ \ cav 

(F±F±) / (F±F±) , will be much smaller than 1 

(^opt± ^ ^opt = (^/41])^). Moreover, we assume 
Tm^th <^ ^7 Tm <^ Topt and ^ <C 1^ as required for 
ground state cooling. The full forces F± = F^^ + F^ 
also contain thermal fluctuations, F^, independent from 
^c^av low mechanical damping (F^ <C 0.) the spec- 
trum of the thermal bath can be replaced by the values 
at the resonances, i.e., {F±F±)l ^ (^±^±)Lsgn(a;)f^± • 
The general scheme yields the variances by integrating 
{v±V±)u = {{F±F^)l + (F±F±)--) I x±H P and 

In the optimal cooling regime, the variances of the 
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7. CONCLUSIONS AND OUTLOOK 



optomechanically damped system can be expressed in a 
compact way: 

2{7r^)/mn± 1 + 2(neff + Sn), (5) 
2mQ±{fil) ^ 2{7Tl)/mQ± + g'^/^l, (6) 

where ngff = F^nth/Fopt + ^opt and Sn = T^nth//^. 

Together with Eq. ([T]), these formulas constitute our 
main result for entanglement in a system subject to op- 
tomechanical cooling. We now extract and discuss its 
main physical features. The first term on the rhs. of 
Eq. ([5| describes the ground state energy, and the sec- 
ond term takes account of the cooling mechanism: the 
thermal occupation is reduced to an effective phonon 
number neff. Thus, entanglement can in principle be cre- 
ated even for large bulk temperatures, nth ^ 1, if the 
optomechanical damping rate Fopt is sufficiently large. 
Since Fopt = 4^^/^, this can be achieved either by re- 
ducing the cavity linewidth tz or by increasing the cavity- 
oscillator coupling rate g. However, we identify two pro- 
cesses that destroy entanglement for small and large 
^, respectively. First, as known from ^16^, the cooling 
mechanism becomes less efficient in the strong coupling 
regime Fopt ^ where the contribution of 6n becomes 
appreciable. Second, for a large optomechanical coupling 
strength ^, the low- frequency contributions of the pho- 
ton shot noise induce an increase of the position variance 
[second term on the rhs. of Eq. ([6|]. This implies that 
strong correlations between the individual oscillators and 
the driven cavity lead to a destruction of entanglement 
between the oscillators. 

As a consequence, entanglement depends non- 
monotonically on the cavity linewidth and the optome- 
chanical damping rate Fopt in the optimal cooling regime, 
cf. Figs. [4] and [5] Entanglement can be generated only if 
the mechanical coupling rate exceeds a value of 

G^in/^ ^ 2(neff + 6n) + Topt^/SQ^ (7) 

Note that Eq. ^ can be employed to optimize entangle- 
ment. 

6. SHORTCOMINGS OF THE LINDBLAD 
APPROACH 

The crucial destruction of entanglement by strong dis- 
sipation is missed entirely by the commonly employed 
Lindblad master equation approach. Its general form 
is given by p = -i[Hsys, p] + ^i^i{p)^ where the in- 
fluence of the bath is taken into account by Lindblad 
terms Ci pjj. For equilibrium baths, these are given 
by = (F^/2)(n± + l)D[i±] and = 

(F^/2)n±P[iy, where V[A]{p) = 2ApA^ - A^ Ap - 
pA^A and A± are the mechanical normal mode anni- 
hilation operators. At zero temperature, the shortcom- 
ings of the Lindblad approach are most obvious: The 




optomechanical damping rate Fopt 



Figure 4. Density plot of the entanglement (logarithmic nega- 
tivity En) for the non-equilibrium dissipative system of Fig.js] 
as a function of the optomechanical cooling rate Fopt and the 
mechanical coupling rate G. The white solid line represents 
the boundary between entangled and separable states of the 
system and thereby defines the minimal coupling rate Cmin 
necessary to observe entanglement. The dashed white line de- 
picts the result for Gmin from the simpler Lindblad approach. 
F^nth 10"^ A± -Q± and k 0.067 

system evolves into its ground state, whose entangle- 
ment is not reduced at all by the system-bath cou- 
pling. To treat the non-equilibrium case of Fig. |3] in the 
Lindblad approach, we have to consider four additional 
terms, = {F±F±y^^l{n±£l/2n)V[A±]{p) and 

45 (p) = (F±F±)!.^^^ (n±£^/21^)D[4](p), which take 
account of the decoherence via the cavity modes (see [19] 
for a detailed derivation). The steady-state variances of 
the normal modes follow as 2m^}±{fj'^) = 2(7r^)/ml]± ^ 
2neff + 1. This expression describes the cooling to an ef- 
fective phonon number neff but fails to capture the loss 
of entanglement for strong optomechanical coupling (see 
the dashed curve in Fig. [4|. The shortcomings of this 
approach can be understood by noting that the Born- 
Markov approximation, which assumes the bath to have 
a very short correlation time (no memory) and to be un- 
correlated with respect to the system, does not hold in 
general for a non-equilibrium bath as can be seen in our 
example. 

7. CONCLUSIONS AND OUTLOOK 

The general exact framework introduced here can be 
employed to analyze the entanglement of oscillators un- 
der the influence of arbitrary bath spectra, among them 
non-equilibrium and tailored non-standard spectral den- 
sities. As pointed out in this paper, the effects of tunable 
non-equilibrium environments promise rich physics to be 
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Figure 5. Minimal coupling strength Cmin necessary to ob- 
serve entanglement. Cmin displays an optimum for interme- 
diate values of the optomechanical damping rate Fopt and the 
cavity decay rate n (for Fm^th = 10~^ 11, A± = — ^±). 



explored in current experimental setups. The optome- 
chanical setup investigated here is in fact just one of a 
rather large class of setups to which this work applies, 
and which also extends into the fields of electromechan- 
ics [5, 6 and cold-atom physics [18 . We also note that 
completely different systems show similar entanglement 
production effects under non-equilibrium conditions, as 
has been explored in the case of coupled, driven qubits 
[20^ , atoms [21^ and ions [22^ , or coupled double quantum 
dots [23]. 

In the quest to observe entanglement in dissipatively 
cooled optomechanical or nanoelectromechanical systems 
the theory presented here serves as an essential guideline: 
It identifies viable parameter regimes for generating and 
optimizing entanglement between massive mechanical os- 
cillators. 

Recent works [191 [24] have proposed an alternative way 
of generating entanglement in nanomechanical systems: 
By modulation the coupling strength between the oscil- 
lators the system can be parametrically driven into a 
non-equilibrium state which features entanglement even 
at relatively large temperatures. In a future work the 
general framework introduced here can be employed to 
discuss the generation of entanglement in a parametri- 
cally driven system and to compare and connect the two 
approaches. 
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